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Abstract 


Differential  equations  are  presented  that  describe  wall 
oscillations  of  a  gas-pressurized  thin-walled  cylinder 
for  cases  where  variations  in  the  axial  direction  are 
negligible.  A  linear  stress/strain  relation  is  assumed. 
Harmonic  solutions  are  obtained.  Limiting  forms  of 
these  solutions  are  given  for  cases  where  the  mode 
number,  n  =  0,  1,  2,...,  is  moderate  and  for  cases  where 
the  mode  number  is  large.  In  the  former  case,  curvature 
effects  are  important.  In  the  latter  case,  the  vibrations 
behave  locally  like  those  on  a  flat  plate.  The 
applicability  of  the  harmonic  solutions,  for  evaluating 
the  hoop  stress  induced  by  a  high  energy  pulsed  laser 
beam,  is  discussed.  The  maximum  stress  induced 
during  the  initial  transient  is  investigated.  An  estimate 
is  given  for  the  maximum  stress  perturbation  induced 
by  a  “slab”  type  (no  axial  variation)  beam  with  uniform 
fluence  and  a  width  equal  to  the  cylinder  diameter.  In 
this  case,  the  estimate  for  the  maximum  value  of  the 
laser  induced  hoop  stress  perturbation,  ,  is  found 
fro™  o^/Oj  =  (CF/pa)(E'/p)‘'^  where,  0;,C,  F,  p, 
a,  E' ,  and  p  are  initial  pressure  induced  hoop  stress. 


coupling  coefficient,  fluence,  cylinder  internal  pressure, 
cylinder  radius,  effective  wall  modulus  of  elasticity, 
and  wall  density,  respectively.  The  maximum  stress 
perturbation  induced  by  slab  and  circular  cross-section 
laser  beams,  with  non-uniform  profiles  and  widths 
which  are  small  compared  with  cylinder  radius,  are  also 
estimated.  For  the  case  of  a  narrow  slab  beam,  the 
upper  bound  on  the  laser  induced  hoop  stress  is  found 
from  am/cii  =  0.5(CFo/pa)^  (E  /p)  where  Fq  is  center 
line  fluence. 


1.  Introduction 


The  interaction  between  a  high-energy  laser  pulse  and  a 
surface  can  result  in  material  blow  off  and  impulse 
loading.  The  latter  may  induce  vibration  and/or 
structural  damage.  As  a  result,  high-energy  pulse  lasers 
are  often  considered  as  a  component  of  a  missile 
defense  system. ' 

Research  in  the  area  of  pulsed-laser  surface  interaction 
has  focused  on  the  determination  of  the  laser  levels 
needed  to  induce  blow-off  and  determination  of  the 
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coupling  coefficient  that  relates  incident  fluence  to 
impulse.  *  The  subsequent  effect  on  the  structure  has 
received  less  attention.  In  this  connection,  Sutton^  has 
evaluated  the  wall  vibration  induced  in  a  pressurized 
cylinder  by  a  high-energy  laser  pulse  of  uniform  width. 
However,  Sutton’s  model  neglects  curvature  effects, 
which,  in  fact,  play  a  significant  role  in  the  structural 
response  for  beam  widths  of  the  order  of  cylinder 
diameter.  Hence,  a  study  of  the  vibrational  response  of 
a  gas  pressurized  cylinder,  to  a  high-energy  laser  pulse, 
including  curvature  effects,  was  undertaken.  A  linear 
stress/ strain  relation  is  assumed.  Emphasis  is  placed  on 
the  case  treated  by  Sutton^,  namely,  an  incident  laser 
beam  of  uniform  fluence  with  width  equal  to  the 
cylinder  diameter  and  infinite  axial  extent.  (Beams  of 
infinite  axial  extent,  with  a  fluence  profile  that  does  not 
vary  with  axial  position,  are  referred  to  as  “slab” 
beams,  herein).  Slab  and  circular  cross-section  laser 
beams,  with  widths  that  are  small  compared  to  the 
cylinder  diameter,  are  also  considered. 

The  response  of  the  cylinder  is  expected  to  consist  of  an 
initial  transient  followed  by  harmonic  motion.  In 
Section  2  (Analysis),  equations  are  presented  which 
describe  the  cylinder  wall  motion  in  the  limit  of 
negligible  end  wall  effects.  Harmonic  solutions  are 
obtained.  In  Section  3  (Applications),  the  applicability 
of  these  harmonic  solutions,  for  evaluation  of  laser 
induced  vibration,  is  discussed.  The  initial  transient, 
caused  by  the  laser  pulse,  is  also  discussed.  Estimates 
are  given  for  the  initial  maximum  hoop  stress 
perturbation  induced  by  the  laser  pulse.  A  comparison 
of  the  present  results  with  those  of  Sutton  is  given  in 
Section  4.  The  effect  of  cylinder  translation  is 
discussed  in  Section  5. 


Analysis 


We  consider  wall  vibrations  of  a  gas  pressurized 
cylinder  of  radius  a,  axial  length  L,  thickness  h,  internal 
pressure  p,  and  wall  density  p  (Fig.  1)  and  assume 
h/a  «  1  and  a/L  «  1 . 


Figure  1 .  Notation  for  evaluation  of  wall  vibrations 
induced  in  pressurized  cylinder  due  to 
pulsed  laser  beam. 

The  initial  pressure-induced  hoop  stress  aj  is 


Let  w  and  v  denote  small  displacements  in  the  radial 
and  tangential  directions,  respectively,  of  the  center  of 
an  arc  element.  These  are  functions  of  the  initial 
angular  position  0  and  time  t,  and  define  the 
displacement  of  the  centroid  of  each  arc  element  from 
its  initial  (t  =  0)  position.  Variations  of  w  and  v  in  the 
axial  direction  are  ignored.  This  is  consistent  with  the 
neglect  of  end  wall  effects  (i.e.,  the  assumption  a/L«l) 
and  the  assumption  of  a  slab  beam  with  a  fluence 
profile  that  does  not  vary  in  the  axial  direction.  The 
neglect  of  axial  derivatives  in  Eqs.  4  of  Fung^  leads  to 
momentum  conservation  equations  given  by 
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Here,  v,  E,  and  E'  =  E/(l-v^)  are  Poisson’s  ratio,  the 
modulus  of  elasticity  and  an  “effective”  modulus  of 
elasticity,  respectively.  The  quantity  ©o  is  the 
fundamental  cylinder  vibration  frequency  and  gj  is  a 
parameter  which  characterizes  the  initial  radial  strain 
induced  by  pressurization.  For  cases  where  the 
separation  distance  between  the  end  walls  is  kept 
constant  during  the  pressurization  process,  the  axial 
strain  is  zero  and  Ei  is  the  actual  radial  strain,  w/a.  If 
the  end  walls  are  unconstrained  during  pressurization, 
an  axial  stress  equal  to  half  the  hoop  stress  is  induced, 
and  the  corresponding  radial  strain  is  then 
Wj/a  =  0.5(2-v)Ei/(l-v^). 


Harmonic  solutions  of  Eqs.  2  are  now  obtained. 
Assume  perturbations  of  the  form 


k.  -  1  +  £.  +  — 
12 


-  1  n 
Va, 


(5a) 


kj  =  1  + 


12 


h 

Vay 


(5b) 


,  14-  2  1  [  h  I  4 

k,  =  1  +  £  n  + —  —  I  n 
12l  a 


(5c) 


1  (  h  1  2 
k4  =1  +  —  —  I  n 

12l  a 


(5d) 


The  kinetic  energy  in  each  mode,  at  t  =  0,  is 

2(ke)„ 


Ttphaffljj 
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v„  ' 


(6) 


Equation  6  defines  the  net  (kinetic  plus  potential) 
energy  associated  with  the  vibrational  motion.  Let  cj 
and  E  denote  the  perturbation  in  the  tensile  (hoop)  stress 
and  in  strain.  These  are  related  by 


w  =  ZWjj  cos  n0  sin  ©^^t 


(3a) 


o  =  E'e 


(7) 


V  =2  VjjSinn0sin©jjt 


An  expression  for  e,  in  terms  of  v  and  w,  is  given  by 
^  '  Eq.  A-2,  namely 


where  n  is  the  vibrational  mode  number  which,  for  each 
mode,  equals  the  number  of  sine  waves  around  the 
cylinder  circumference.  Substitution  of  Eqs.  3  into  Eqs. 
2  yields 
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Further  discussion  relating  to  the  computation  of  strain 
is  given  in  Appendix  A. 

Simplified  equations  for  the  frequency,  velocity  ratio, 
and  strain  associated  with  the  positive  and  negative 
roots  of  Eq.  4a  are  now  given.  We  make  the  realistic 
assumption 
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If  terms  which  are  small,  compared  with  one,  are 
neglected,  the  dependent  variables  associated  with  the 
positive  root  of  Eq.  4a  become 

(®n/®o)  =  1  +  n^  (10a) 

Vn/wn=n  (10b) 

=  (l  +  n^)^wjj/a  j  cos  n0  sin  coj^t  + 


These  are  the  equations  that  describe  the  fundamental 
vibration  mode  for  a  flat  plate  of  width  s  =  tia/n.  This  is 
a  consequence  of  the  fact  that,  for  large  n,  the  arc  length 
between  radial  velocity  nodes,  namely  7ta/n,  is  small 
compared  to  the  cylinder  radius  a  and  thus  departs  only 
slightly  from  a  flat  plate. 

Equations  10c,  1  Ic,  and  12c  define  the  strain  associated 
with  a  single  mode.  In  the  case  of  multiple  modes,  a 
summation  process  is  needed  as  discussed  in  Appendix 
A. 
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Similarly,  the  dependent  variables  associated  with  the 
negative  root  of  Eq  4a  become 


(11a) 


Vn/wn  =  -1/n 


(11b) 
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Note  that  cOj^  =  0  for  n  =  0,l. 

It  can  now  be  shown  that  the  positive  and  negative  roots 
of  Eq.4a  have  the  following  physical  significance. 
Substitution  of  Eq.  10b  into  Eq.  2b  indicates  that  is 
a  consequence  of  the  conservation  of  v  (transverse) 
momentum.  Similarly,  substitution  of  Eq.  1  lb  into  Eq. 
2a  indicates  that  is  a  consequence  of  the 
conservation  of  w  (radial)  momentum. 


In  the  limit  of  large  n,  namely  1/n  «  1,  Eqs.  1 1  become 


3.  Applications 

We  now  investigate  the  stress  perturbation  induced  by 
an  energetic  laser  pulse  of  width  s,  incident  on  the 
cylinder  surface.  It  is  assumed  that  ablation  creates  an 
impulsive  load  which  results  in  inward  momentum  of 
the  cylinder  surface.  The  resulting  stress  perturbation  is 
estimated. 

Consider  the  case  of  a  laser  pulse  of  uniform  fluence  F 
and  width  s  (Fig.  1)  incident  on  the  cylinder  surface  at 
time  t  =  0"^.  Let  C  represent  the  coupling  coefficient, 
namely  the  ratio  of  the  induced  impulse  (radial 
momentum/area)  to  the  incident  fluence  (energy/area). 
The  initial  conditions  regarding  cylinder  surface  motion 
are  then,  for  t  =  0^, 
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where  0s  is  the  value  of  0  corresponding  to  the  edge  of 
the  laser  beam, 

0g  =  sin  ’ [s/(2a)]  (13d) 
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12  ' 
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The  initial  condition  v  =  5v/0t  =  0  is  clearly 
inconsistent  with  Eqs.  3  and  4b  so  that  those  equations 
cannot  be  used  directly  to  construct  a  solution  which 
satisfies  the  initial  condition  v  =  0,  except  for  the  case 
of  large  n.  (See  Eq.  12b).  There  will  be  an  initial 
transient  in  the  region  |0|  <  9^  which  later  spreads  to 
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encompass  the  entire  cylinder  and  evolves  into  a 
multimode  harmonic  motion.  It  is  expected  that  the 
peak  induced  stress  will  occur  during  the  initial 
transient  [times  of  order  cot  =  0(1)].  An  estimate  of  this 
peak  stress  is  given  below  for  the  cases  s  =  2a  and  s/a 
«  1. 

Case  s  =  2a: 


In  view  of  the  initial  condition  v  =  0,  an  estimate  of  the 
peak  stress,  during  the  initial  transient,  can  be  found 
from  Eq.  2a  by  assuming  v  =  0,  therein,  and  by 
assuming  that  the  vibration  remains  confined  to  the 
region  |0|  <  nj 2.  Thus  we  assume  cylinder  wall 
perturbations  of  the  form 


cosOsincot 

(14a) 

phco 

v  =  0 

(14b) 

=  w/a 

(14c) 

for  |0|  <  71/2.  Substitution  into  Eq.  2a,  and  neglect  of 
terms  of  order  sj  and  (h/a)^/12  yields 

2  2 

CO  =coq  (15) 

The  maximum  strain  gjjj  occurs  at  0  =  0  and  cot  =  37i/2, 
and  equals 


CF 

ph 


vE'y 
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the  quantities  Wj^  and  co,  .  The  superscript  plus  sign 
is  henceforth  omitted.  The  initial  kinetic  energy  is 
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The  energy  in  mode  n  =  1  is  (Eq.  6) 
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Equating  Eqs.  17  and  18,  and  recalling  Eq.  10,  indicates 


Wj 

a 


C¥f  1 


ph 


E' 


(19a) 


2 


COj 


(19b) 


The  maximum  strain  and  stress  are  given  by 
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Equation  20b  is  only  a  factor  1/ yjl  less  than  the  value 
in  Eq.  16b.  The  degree  to  which  Eq.  20b  characterizes 
the  late  time  oscillation  requires  further  study. 


Case  s/a  « 1: 


which  may  be  viewed  as  providing  an  upper  limit  on 
the  induced  strain.  The  maximum  hoop  stress,  due  to 
the  laser  pulse,  is  found  from  =  E'  8^^^,  or 


CF 

pa 
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1/2 


V  p  y 
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An  estimate  of  the  stress  perturbation  associated  with 
the  long  term  harmonic  motion  is  now  found.  We 
assume  that  Eqs.  3  to  10  apply  and  that  the  induced 
vibration  is  primarily  the  n  =  1  mode  with  energy  equal 
to  the  kinetic  energy  deposited  at  t  =  0.  This  choice 
most  closely  matches  the  initial  radial  velocity  profile 
for  the  region  |0|  <  njl,  and  corresponds  to  the  use  of 


In  the  present  case,  the  laser  spot  size  is  small,  relative 
to  the  cylinder  radius,  and  a  flat  plate  approximation 
can  be  made  (Fig.  2).  The  transverse  coordinate  v  is 
now  replaced  by  the  Cartesian  coordinate  y,  as 
indicated  in  Fig.  2.  It  is  assumed  that  the  perturbation 
is  initially  confined  to  the  spot  area  (|y|  <  s/2)  and  we 
obtain  an  estimate  of  the  corresponding  maximum 
strain.  This  estimate  is  an  upper  limit  since  the 
disturbance  spreads  laterally  and  thereby  attenuates. 
We  consider  both  a  “slab”  beam  (beam  fluence  profile 
independent  of  axial  position)  and  a  beam  of  circular 
cross-section.  The  late  time  harmonic  motion 
associated  with  a  narrow  slab  beam  is  discussed  in 
Appendix  B. 
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Oj  ^  ®i 

-s/2  s/2 

Figure  2.  Notation  for  case  of  slab  laser  beam,  or 
circular  cross-section  beam,  incident  on 
slightly  curved  cylinder  section  (d/w  «  1). 


Slab  Beam:  Using  a  flat  plate  approximation,  the 
fundamental  frequency  and  mode  shape  is  of  the  form 


s/h  be  of  the  order  10,  or  more.)  The  induced  vibration 
is  then 
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The  maximum  induced  strain  occurs  at  y  =  s/2  and 
equals 
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V  9y  y 


1  (CTq)' 

2  ph  pa 


(25a) 


For  a  given  cylinder  and  fluence  Fq,  the  induced 
maximum  strain  is  independent  of  slab  beam  width  and 
is  therefore  independent  of  laser  beam  energy.  The 
corresponding  induced  stress  is 
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w~  cos(7iy/s)sin<ot  (21b) 

Equations  21a  and  21b  are  obtained  by  introducing 

C0  =  C0q,  W  =  Wjj,  s  =  7ia/n and y  =  0a  into  Eqs.  3a 

and  12.  For  convenience,  consider  an  incident  fluence 
profile  of  fhe  form 


F 


(22a) 


with  energy  per  unit  axial  length  equal  to 

T 

L  -s/2 


Y  s/2 
-=  I  Fdy= 


2F„s 


(22b) 


Flere,  J  denotes  net  beam  energy.  The  vertical  velocity 
at  time  t  =  0  is 


5w 

at 


CFq 

(-1)  — ^cos 
ph 


(23) 


(Note  that  Eq.  23  does  not  apply  for  values  of  2y/s  near 
one,  since  the  magnitude  of  the  fluence  in  the  region  is 
below  the  threshold  needed  to  induce  ablation. 
Flowever,  this  region  has  a  small  impact  on  the  induced 
strain.  The  use  of  the  product  CF  to  characterize  the 
local  impulse  given  to  the  wall  probably  requires  that 


Equations  25a  and  25b  may  be  viewed  as  providing 
upper  bounds.  Lateral  expansion  reduces  the  vibration 
amplitude.  An  estimate  of  this  effect  is  found  as 
follows.  For  the  case  of  large  n  and  negligible 
thickness  effects,  Eq.  10  has  a  solution  of  the  form  w  ~ 
f(a9  -  ct)  -I-  g(a0  -I-  ct)  which  represents  waves  with 
velocity  c  =  (a,  /p)'  ^  traveling  in  the  -1-0  and  -0 
directions,  respectively.  After  a  quarter  period,  cot  = 
n/2,  the  disturbed  area,  originally  of  width  s,  has 
expanded  to  a  width  of  2s.  Equation  25a  indicates  that, 
for  fixed  energy  (i.e.,  sFq  =  consfant),  the  maximum 
strain  is  inversely  proportional  to  the  square  of  the 
beam  width.  Thus,  the  maximum  strain  at  cot  =  nil  can 
be  estimated  from 


_  1  (CFpf 
8  phpa 


(25c) 


which  is  !4  the  value  in  Eq.  25a.  The  strain,  and 
corresponding  stress,  decrease  with  further  expansion 
of  the  vibrating  region. 

Circular  cross-section  beam:  We  now 
consider  a  circular  cross-section  beam  with  radius  s/2. 
Here,  the  quantity  y  in  Fig.  2  corresponds  to  a  radial 
coordinate.  In  order  to  estimate  the  initial  induced 
stress,  in  the  present  case  of  cylindrical  symmetry,  it  is 
necessary  to  consider  the  pressure -induced  axial  stress, 
as  well  as  the  tangential  stress  ai.  For  unrestrained  end 
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walls,  the  pressure-induced  axial  stress  equals  qJ2. 
Consistent  with  the  objective  of  estimating  an  upper 
limit  on  the  induced  stress,  we  will  assume  that  both  the 
axial  and  the  transverse  pressure -induced  stresses  equal 
ai  and  neglect  thickness  effects.  In  this  case  the 
induced  vertical  perturbations  are  axisymmetric  and  are 
described  by  the  “membrane”  vibration  equation 


,2 

p  5  w 


1  d  f  0w  ^ 

y  3y  V  dy  J 


(26) 


which,  for  zero  displacement  at  y  =  s/2,  has  a  solution 
of  the  form 


2.405 


sincot 


CO  =  2.405 


^a.V/2 


V  p  y 


(27a) 


(27b) 


with  limitations  similar  to  those  observed  for  Eq.  23. 
The  corresponding  vibration  and  induced  strain  are 


,  .  CF„  (  2y  , 

(-l)  - Jg  2.405 —  I  sincot  (31a) 

phco  V  s 


2  \  5y  /  2 


ph 


0  I  P  2y 


—  Jj  I  2.405  — 


(31b) 


The  maximum  strain  occurs  at  2.405  (2y/s)  =  1.841  and 
equals 


8 


m 


0.1693 


(CFq)' 

ph  pa 


(32a) 


Here  again,  the  maximum  strain  depends  on  Fg  and  is 
independent  of  beam  radius  and  beam  energy.  The 
corresponding  stress  equals 


where  Jm(x)  is  o  Bessel  Function  of  the  first  kind,  of 
order  m.  We  note 


0.1693 


L  pa  J  p 


(32b) 


Jo  (x)  =  (-l)Ji  (x)  ,  Jo  (2.405)  =  0  (28a) 

Jj(l.84l)  =  0.5819(  max  pt)  (28b) 

2.405 

j  xJg(x)dx=  1.2485  (28c) 

0 

For  convenience,  we  consider  an  incident  beam  of  the 
form 


Equations  32a  and  32b  represent  upper  bounds.  Lateral 
expansion  decreases  the  strain.  An  estimate  of  this 
effect,  for  circular  beams,  is  obtained  as  follows. 
Assume  that  the  edge  of  the  disturbance  moves  with  a 
radial  velocity  equal  to  c=  (oj/p)'  .  After  a  quarter 
period,  the  beam  diameter  increases  from  s  to  s(l  -l- 
71/4.810).  Equation  32a  indicates  that  for  fixed  energy 
(i.e.,  s^Fq  =  constant),  the  maximum  strain  varies 
inversely  as  the  fourth  power  of  the  diameter.  Hence 
the  maximum  strain  at  ot  =  n/2  can  be  estimated  from 


—  =  Jq  [  2.405  —  I  (29a) 

Fg  I  s  J 


0.1693  (CFq)^ 
7.469  ph  pa 


(32c) 


which  contains  a  net  energy  of 
S/2 

T  =  2;i  j  Fydy  =  0.4317  ;i 
0 


The  induced  vertical  velocity  at  t  =  0  is 


0w  ,  .  CF  /  2y 
—  =  (-l)  —  Jo  2.405  — 
5t  ph  V  s 


which  is  a  factor  1/7.469  less  than  the  value  in  Eq.  32a. 
The  strain  decreases  with  further  expansion  of  the 
vibration  area. 

(29b) 

Equations  25b  and  32b  indicate  that,  for  fixed  values  of 
the  independent  variables,  the  maximum  stress  induced 
by  a  narrow  slab  beam  is  a  factor  2.95  larger  than  that 
induced  by  a  corresponding  small  diameter  circular 
cross-section  beam.  This  suggests  that  the  use  of  the 
(30)  narrow  slab  beam  solution  to  estimate  the  strain 

induced  by  a  small  diameter  circular  cross-section  beam 
(as  is  done  in  Ref  4)  will  lead  to  an  over  estimate  of 
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induced  strain  by,  approximately,  a  factor  2.95.  The 
latter  is  approximate  in  view  of  the  assumption  that  the 
initial  axial  stress  equals  the  initial  tangential  stress,  Gi, 
in  the  small  diameter  circular  cross-section  beam  case. 


4.  Comparison  with  Sutton 

Sutton^  considers  a  uniform  incident  beam  with  a  width 
of  the  order  of  the  cylinder  diameter.  He,  in  effect, 
applies  flat  plate  equations  (i.e.,  the  equivalent  of  using 
Eqs.  12)  and  thereby  neglects  curvature  effects.  The 
resulting  solution  can  be  deduced  by  taking  s  =  2a  in 
Appendix  B.  Sutton  excludes  the  n  =  1  mode  and 
claims  that  the  n  =  2  mode  is  the  major  contributor  to 
cylinder  strain.  His  maximum  strain  results,  in  present 
notation,  are  (see  Eq.  B-8) 


^  CF^ 


V37iy 


ph  pa 


(33a) 


=  0.02252 


VPaj 


E' 

P 


(33b) 


We  now  compare  the  fluence  required  to  produce  a 
given  value  of  Gm/ci  as  given  by  the  present  model 
(Eq.  16b)  with  the  value  obtained  from  the  Sutton 
model  (Eq.  33b).  Let  Fpj-gsent  denote  the  fluence  in 
Eq.  16b  and  let  Fgutton  denote  the  fluence  in  Eq.  33b 
The  ratio  of  these  quantities  is 


of  the  quadratic  terms  in  Eq.  8  yields  reduced  estimates 
of  strain  and  incorrect  scaling  laws. 

In  a  more  recent  study"*,  Sutton  repeats  the  above 
solution  for  the  case  s  =  2a.  He  also  includes  a  study  of 
the  case  s  «  2a,  for  which  the  flat  plate  assumption  is 
valid.  Sutton  claims  that,  for  the  case  s  «  2a,  the  n  =  2 
mode  is  the  dominant  mode  with  respect  to  induced 
strain.  This  claim  appears  to  be  incorrect  for  reasons 
discussed  in  Appendix  B. 

5.  Cylinder  Translation 

The  cylinder  is  assumed,  initially,  to  be  stationary  in 
space.  The  laser  pulse  at  t  =  0  will  induce  a  uniform 
translation  of  the  cylinder  center  of  mass.  The  solution 
of  the  present  problem,  in  a  coordinate  system  which 
moves  with  the  center  of  mass,  is  discussed  herein. 

Consideration  of  conservation  of  linear  momentum 
indicates  that  the  laser  pulse  will  induce  a  cylinder 
center  of  mass  velocity  V,  in  the  0  =  180°  direction, 
equal  to 


ph  20„  -I-  sin20, 

—  ^  (35a) 

CF  471 


=  1/4  0s  =  7t/2  (35b) 


=  0s/7t  0s«l  (35c) 


p 

sutton 

V 

present 


(  1 
^0.02252  Gjjj/gj  ^ 


(34a) 


Note  that,  for  0s  =  n/2,  the  magnitude  of  V  is  !4  the 
value  of  dw/dt,  at  0  =  0,  in  Eq.  13b.  Also,  V  becomes 
vanishingly  small  for  small  width  slab  beams  (0s«l). 


=  21.07  G^/Gi  =  0.1  (34b) 

=  12.17  a^/G-0.3  (34c) 

Equations  33  overestimate,  by  an  order  of  magnitude, 
the  fluence  required  to  achieve  a  maximum  stress  of 
order  Gm/cri  =  0  (0. 1).  It  should  be  noted  that  Eq.  16 
refers  to  the  initial  transient  whereas  Eq.  33a  refers  to 
the  subsequent  harmonic  motion.  However,  a 
comparison  of  Eq.  33a  with  Eq.  20b  would  introduce  a 
factor  of  only  1/  x/l  on  the  right-hand  side  of  Eq.  34 
and  thereby  would  not  significantly  affect  the 
conclusion.  In  summary,  the  major  deficiencies  of  Ref 
23,  from  the  viewpoint  of  the  present  study,  are  (a) 
curvature  effects  are  neglected  and  (b)  the  exclusive  use 


Equations  2  are  applicable  in  an  inertial  frame.  The 
initial  conditions  for  Eqs.  2,  in  a  coordinate  system 
fixed  with  respect  to  the  moving  cylinder  center  of 
mass,  become 


II 

<; 

II 

o 

(36a) 

0w  f  cv'] 

—  =  V  -  —  COS0 

a  V  phy 

|0| 

<0S 

(36b) 

=  V  COS0 

0S< 

0|<7I 

(36c) 

cV 

—  =  (-l)V  sin0 

l0l 

<71 

(36d) 
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Unlike  Eqs.  13,  Eqs.  36  correspond  to  zero  net  linear 
momentum.  The  solution  of  Eqs.  2,  with  Eqs.  36  as 
initial  conditions,  defines  the  perturbations  with  respect 
to  the  moving  cylinder. 

Equations  2,  with  Eqs.  13  as  initial  conditions,  are 
convenient  for  finding  the  perturbations  (from  the 
initial  cylinder  location)  at  small  times.  The  condition 
that  w  and  v  remain  small  becomes  violated  with 
increase  in  time.  However,  the  solution  of  Eqs.  2,  with 
Eqs.  36  as  boundary  conditions,  is  valid  for  all  times, 
t  >  0,  and  is  convenient  for  finding  the  late  time 
solution.  The  latter  is  not  explored  in  the  present  study. 

6.  Concluding  Remarks 

The  present  study  provides  an  upper  bound  to  the 
maximum  hoop  stress  induced  in  a  thin-walled  elastic- 
pressurized  cylinder  by  a  high-energy  laser  pulse. 
Further  effort  is  needed  to  refine  these  estimates  and  to 
establish  the  longer  term  vibrational  behavior  of  the 
cylinder.  In  addition,  end-wall  effects  may  need  to  be 
considered  in  some  cases.  Finally,  the  relationship 
between  the  present  linear  stress/sfrain  estimates  and 
structural  failure  mechanisms,  in  pressurized  cylinders, 
needs  further  study. 
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Appendix  A:  Strain  Expressions 

Strain  is  related  to  the  perturbations  v,  w  in  the 
following  section.  Let  aA0  denote  an  elementary  arc 
section. 


a(1  +  £)A0 


Figure  A-1.  Notation  used  for  evaluation  of  strain  s. 

The  strain  induced  by  a  displacement  can  be  expressed 
(Figure  A-1) 


[(l  -l-s)aA0]^  =  (Aw)^ 


,  ,  5v 

(a  +  w)A0  H - A0 

00  . 


n2 


(A-1) 


Expansion,  retention  of  leading  terms  and  application  of 
the  limit  A^  — >  0  indicates 
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In  the  case  of  harmonic  motion,  the  linear  terms  in  Eq. 
A-2  dominate  for  moderate  values  of  n  and  the 
quadratic  terms  dominate  for  large  values  of  n. 


A  summation  is  required  when  there  are  multiple 
modes.  The  strain  associated  with  moderate  values  of  n 
is  (recalling  that  nVj^/wj^  =  -1  for  this  case) 

£  =  z(l  +  n^j  (^wJi/aj  cos  nO  sin  cOj^  t  (A-3) 
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The  maximum  strain  occurs  at  0  =  0.  For  large  n, 
where  the  fiat  plate  approximation  applies,  the  strain  is 
found  from 

-|2 
nw 

Z - ^  sin  nO  sin  t  (A-4) 

n  a 
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If  there  is  a  dominant  mode  number  n,  the  maximum  The  initial  conditions,  associated  with  the  laser  pulse, 
strain  occurs  at  n0  =  7t/2.  can  be  expressed 


Equation  A-2  defines  the  tensile  strain  perturbation 
associated  with  the  mid  section  of  the  cylinder  wall. 

The  total  strain  at  the  mid  section  is  found  from  the  sum 
S  +  Sj.  The  strain  at  the  outer  surface  of  the  cMinder 
is  found  by  adding  the  bending  strain  Ej,  =  h/yuC. 

Appendix  B;  Fourier  Analysis 

For  the  case  of  an  incident  laser  beam  with  small  width, 
the  flat  plate  approximation  (Eq.  12)  is  applicable.  The 
harmonic  motion,  which  follows  the  initial  transient, 
can  then  be  found  by  a  Fourier  analysis.  This  solution 
is  presented  herein.  Superscript  minus  signs  are 
omitted. 

Let  0g  denote  the  value  of  0  corresponding  to  the  edge 
of  the  incident  beam,  namely, 

0g  =sin'^(s/2a)  (B-1) 


Define  a  function  f(0)  such  that 


f(0)  =  cos  0  1  0 

1  <0, 

(B-2a) 

=  0 

0s 

<  lo 

<  n 

(B-2b) 

The  function  f(0),  expressed  as  a  Fourier  series. 

is 

f(0)  = 

00 

1  COS 

nd 

(B-3a) 

2  n  =  1 

where 

1 

r  71  /  \ 

J^^f(0)  cos 

n0  d0 

(B-3b) 

71 

1 

r  0o 

71 

J  0  COS0  cos  n9  d9 

(B-3c) 

1 

sin  20q 

n  =  l 

=  — 

es  + 

(B-3d) 

71 

L  2  J 

1  f  sin(n-l)0s 

1  sin(n+l)0s^ 

n;^l 

(B-3e) 

71  y  n-1 

n  +  l 

J 

z^dwA  °°  CF 

V  dt  yt=0  n=0  ph 

It  follows  that 


(-1)—  cOjjWjj  =  —  n  =  0 
CF  2 


(B-5a) 


-  cL^  n  =  l,2, ...  (B-5b) 


The  flat  plate  approximation  for  strain  is 


E  = 


2a 


2]  nWjj  sinnOsincOjjt 
.n=2 


(B-6) 


Recall 


£:CO„ 


p/(pha) 


and,  for  »1  and  ^n^/l2 j(h/a)^  «  1, 
(On  =  ncooJ^ 


(B-7a) 


(B-7b) 


The  strain  is  then  given  by 


2phpa 

{erf 


2]  a^j  sin  n0  sin  co^jt 
.n=2 


(B-8) 


Equation  33a,  in  the  body  of  the  report,  can  be  obtained 
by  taking  n  =  2,  0^  =  n/2,  02  =  7i/4  and  (02t  =  n/2  in  Eq. 
B-8.  We  now  assume  that  the  beam  has  a  narrow  width 
so  that 


0s  =  s/(2a)  «1  (B-9) 


It  follows  that  for  moderate  n, 

ajj  =  20g/7i  (B-10) 

Let  N  donate  the  mode  number  for  which  the  arc  length 
between  radial  velocity  nodes  is  equal  to  s,  namely 

N  =  7ia/s  =  7i/(20j  (B-11) 
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For  values  of  n  such  that  n/N=0(l), 


an  =  20J71  (B-12) 


With  further  increase  in  n,  the  magnitude  of  Un 
approaches  zero.  Hence  we  take 

ajj  =  20g/7i  n<N  (B-13) 

Substitution  into  Eq.  B-8  yields  the  following 
expressions  for  strain 


71 

2 

3phpa 

r  N  1 

— 

?  £ 

2^  sm  nO  sinco^t 

UeJ 

(cf) 

_n^2 

(B-14a) 


<  N 


2 


(B-14b) 


The  inequality  in  Eq.  B-14b  follows  from  the 
observation  that  each  term  in  the  summation,  on  the 
right  hand  side  of  Eq.  B-14a,  is  equal  to  or  less  than 
one.  The  maximum  strain  then  satisfies  the  inequality 


(CT)" 

2phpa 


(B-15) 


The  magnitude  and  location  of  the  peak  strain  can  be 
found  from  a  numerical  evaluation  of  the  summation  in 
Eq.  B-14.  Equation  B-15  is  consistent  with  Eq.  25a 


For  moderate  values  ofn,  the  above  solution  does  not 
satisfy  the  initial  condition  v  =  0,  as  can  be  seen  from 
Eq.  1  lb.  However,  in  the  limit  of  small  s/(2a),  these 
modes  have  negligible  impact  on  the  solution.  The 
above  solution  represents  an  asymptotic  limit  following 
the  initial  transient.  The  degree  to  which  this  solution 
is  achieved  depends  on  the  rate  of  dissipation  of  the 
vibrational  modes. 

Sutton"^  uses  a  Fourier  expansion  to  describe  the 
oscillations  induced  by  a  laser  beam  for  both  the  s/(2a) 

=  1  and  the  s/(2a)  «  1  cases.  The  deficiency  of  his 
s/(2a)  =  1  solution  was  discussed  in  the  body  of  the 
report  (Eq.  34).  His  treatment  of  the  s/(2/a)  «  1  case  is 
now  briefly  noted.  In  this  case,  as  in  the  previous  case, 
he  claims  that  the  major  contribution  to  the  strain 
comes  from  the  n=2  mode.  The  corresponding 
maximum  strain,  using  the  present  model,  is 

2  (CF)^  2 

7t  phpa 

Equation  B-16  differs  from  Eq.  25a  by  a  factor 

1^2  0  due  to  the  fact  that,  after  the  initial 

transient,  the  disturbance  is  distributed  throughout  the 
cylinder.  Moreover,  the  Sutton  assumption  that  the  n  = 
2  mode  is  the  dominant  mode,  with  respect  to  induced 
strain,  appears  to  be  incorrect  in  view  of  Eq.  B-14a. 
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